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of the dipole-dipole interaction. 
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1. Introduction 



For fluorescing quantum systems one of the most important statistical quantities is 
given by the intensity correlation function, g(r), for photon counts Q. Its behavior for 
small times r can indicate more classical or more quantum behavior and bunching or 
anti-bunching, depending on whether g(0) > 1 or g(0) < 1. Early investigations of the 
intensity correlation function of single two level systems 0, |3], || led to the observation 
of nonclassical light ]5], ||. The intensity correlation function can also exhibit Rabi 
oscillations, and for V systems with metastable state it contains indications of light and 
dark periods ("blinking") in the fluorescence [[?]]. 

The importance of correlation functions partially stems from the relative ease with 
which they can be experimentally determined, partially due to the fact that the efficiency 
of the photon detector does not enter. Correlation functions were determined for single 
ions in a trap in experiments which were mainly carried out for the detection of quantum 
jumps in the fluorescence || [TU], [TT| 



In the last decade the same was achieved 



for single fluorescent molecules |12], [13], [nj. A theoretical determination of g{r) was 
obtained in closed form for three level systems |T5], IB, [T7J and g(0) was calculated for 
two dipole-dipole interacting two-level systems [|TH|, [B|, |2D| , but for more complicated 
multi-level systems it is usually done by the numerical solution of Bloch equations 
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0, ^T], [22|, |23], Because the intensity correlation function contains information on 
the system parameters, e.g. atomic constants like Einstein coefficients, these parameters 
could be obtained in principle from an experimentally determined g{r) by fitting 
numerically calculated curves to the data. However, in general such a fit procedure needs 
many numerical runs and tends to be highly sensitive to experimental and numerical 
errors. With an algebraic expression a fit would be much easier and more reliable. 
Moreover, if algebraic expression were known for the intensity correlation function for 
multi-level systems one could study its behavior without recourse to many numerical 
runs for different parameter values. In the present paper such algebraic expressions of 
the correlation function will be given for fluorescing quantum systems with light and 
dark periods. As a generalization of a result for V systems in it is deduced from 
this algebraic expression that the correlation function for "blinking" systems with dark 
periods shows a "hump" for values of r larger than the correlation times of the individual 
intensity periods. It is also pointed out that the existence of extended dark periods may 
considerably enhance the amplitude of Bloch oscillations. 

As an application we study two dipole-dipole interacting V systems and determine 
an algebraic expression for the corresponding intensity correlation function g{r). This 
expression is used to study the behavior of g{r) with respect to the strength of the 
dipole-dipole interaction. It is shown that the dipole-dipole interaction has an effect 
for small r which is particularly pronounced for small values of the strong driving. We 
suggest to use this effect to experimentally verify the dipole-dipole interaction for atomic 
distances of a few wavelength of the strong transitions. 



2. Algebraic expressions for intensity correlation functions 

The temporal intensity correlation function for photon counts is defined as follows. In 
the steady state, let G(r) be the joint probability density for detecting a photon both 
at time and r, and let iss denote the steady-state intensity (counts per unit time). 
Then g{r) is defined as [^5] 

g(r) = G(r)// S 2 s . (1) 

We consider now a fluorescing system with n periods of different intensity Ii and mean 
duration Tj, i — 0, . . . , n — 1. A particular period usually corresponds to transitions 
within a simpler subsystem which is easier to treat than the complete system. In the 
following we therefore assume that within a given period i the corresponding intensity 
correlation function gi(r) is known. 

The intensity correlation function g(r) of the larger system can then be determined 
as follows. Let Pi be the probability for the occurrence of period i and let -Pjj (t) be the 
probability to have period j at time r provided one had period i at r = 0. Now, if at 
time the fluorescence is in period i, then, for r CTj, one still is in period i with high 
probability. Therefore, in this case, the joint probability density for detecting a photon 
at both times, and r, is J2i PiGi(r) where Gi(r) is the joint probability density for 



The intensity correlation function of "blinking" quantum systems 



3 



period i. Since for small r one has Ph(t) — 1 this can be replaced by 

E^iMG.W, r«T l . (2) 

i 

For larger r one may end up in a period different from that at time 0. If r is larger than 
the correlation times associated with gi(r) then, for initial period i and final period j, 
the joint detection probability is IjP^(r)Jj. Since for such r's one has Gj(r) = J| one 
can write the complete joint detection probability density as 

G(r)=Y,P i m j (r)G J (r)/I 3 ■ (3) 

ij 

This includes equation @ as a special case since for r small one has Pijij) = for 

In this general formula one still has to express Pi and Pij{r) in terms of the transition 
rates Pij from period i to period j. These transition rates themselves can be expressed 
by Einstein coefficients and Rabi frequencies. One has, in particular, T, = 1/ J^k^iPik- 
The Pijij) are easily seen to obey rate equations, e.g. 

Pn(r) = (-Y,Pik)Pii(r) + P2iPi2(r) + ■ ■ ■ + p nl P ln (r) . (4) 



In general, with the matrix B = (P>ij), 

Bij = pij - 8ij Pik , (5) 



k 

and the matrix 



P(r) = (Ph(t)) (6) 

one has 

P = PB , (7) 
with the initial condition Py(0) = 8y, or 

P(0) = 1 . (8) 
The solution of equation (|7|) with this initial condition can be written as 

P(r) = e B ^ . (9) 



If yU , • • • , A*n-i are the eigenvalues of B (assumed distinct) then 126 



n-l 



B - n c 



3 Br = JJ . (10) 

i=0 a^i Z 1 * 



The properties of the matrix B are closely related to those of stochastic matrices [^g, [27 



and under quite general conditions B has a single eigenvalue /j,q = and eigenvalues 
Hi, . . . , Hn-i with negative real part. This will be assumed in the following. For up to 
four intensity periods, the /ij can be determined in closed form. 

To find the Pj's, we note that for r — » oo the memory to the initial start is in 
general lost. Therefore, for any k, 

P = P K Aoo). (11) 
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(12) 



Laplace transforming equation (0) yields 
P(oo) = lim e(e - B)" 1 . 

Thus the p's can be calculated without knowledge of the //j's. 

The steady-state fluorescence intensity is iss = YL% Pih ■ Therefore we obtain 

9 <r) = ^ = . (13) 

For r larger than the correlation time of gi{r) one has g%{r) ~ 1, and thus equation (O) 
becomes 



(r) = ZijPOiljPi^r) (14) 

fr pi) 

so that for larger r the dependence on r of g(r) is solely governed by the statistics of 
the individual periods. 



strong laser field, 
G 3 




weak field, 
£l 2 



Figure 1. V system with metastable level 2 and Einstein coefficient A3 for level 3. 
^2 and ^3 are the Rabi frequencies of the two lasers driving the weak 1-2 transition 
and the strong 1-3 transition, respectively. 



A first simple example is a system with dark (Iq = 0) and light (I\ 7^ 0) periods, such 
as the V system of figure [1] with a metastable level or a A system []TTJ, |28], |9|, [30], |31], [32 



The mean period durations are T = 1/poi an d T\ = 1/pio, respectively. Because I = 0, 
equation ([13|) becomes 

<?(t) = ^h(t)<?i(t) , (15) 

with gi (r) the correlation function of the two-level subsystem. A simple calculation 
yields 

Pn(r) = ~rp~~T~^p~ + Tf^e-^^ (16) 
J o + J 1 J o + J 1 

J o + 1 1 

For a system equation (p~5|) agrees with the result of [17| if the correlation time for 
gi(r) is much smaller than T\. 

For r <C Ti one has Ph(t) = 1 and g(r) = g\(r)/P\. Since Pi < 1 it follows that, 
for small r, g{r) is just gi(r) blown up by the factor 1/Pi. In particular, possible Bloch 
oscillations of gx{r) become enhanced in g(r) if T increases. Moreover, for r values 
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larger than the correlation time of the two-level subsystem one has g\{r) = 1 so that 
there is a "hump" larger than 1 in g{r). In the case of a V system this hump was already 
noted in 0. 

Such a hump is a general feature for any system with a dark period, as will now be 
shown by means of equation (fH]). Indeed, for r larger than the correlation times of gi(r) 
but much less than the Tj's one has Py = Sy, resulting in g{r) = J2% Pilf/iHa Pal a ) 2 - 
By Schwarz's inequality one obtains (£ a Pi /2 P Q 1/2 / a ) 2 < (E^o P*)(£ a PJl)- Since 
J2 a ^o Pa = 1 — Pq < 1 t ne statement follows. This does not mean, however, that g(r) 
always stays above 1 for all subsequent values of r. There are systems where it dips 
below 1 again and then approaches its asymptotic value 1 from below. 

In the following, equation (|T3"D for g(r) will be studied for an example with three 
different fluorescence periods and shown to be highly accurate. 

3. Application to two dipole-dipole interacting atoms 

We consider two dipole-dipole interacting V systems as in figure ^ at a fixed distance 
r, with one laser driving the strong 1-3 transition and another the weak 1-2 transition. 
Such a system exhibits three fluorescence periods, a dark period Io = 0, a period I\ 
and a double intensity period I 2 — 2I\. The transition rates py between the periods 
are known [Q and are given in the Appendix. One has T = 1/poi; T\ — V(Pio +V12), 
T 2 = l/p2i> and P02 = P20 = 0. The distance dependent complex dipole-dipole coupling 
constant C3 is also given in the Appendix. For r — >• 00 the system behaves as two 
independent, non-interacting, fluorescing V systems, and then the periods I2, I\ and Io 
correspond to two, one or no atom radiating. Because of I = 0, equation flT3"| ) simplifies 
to 

g(r) = {PJlP^g^r) + P 1 hhPi2{r)g 2 {r) + P-zhhP^g^r) + P2llP22{r)g 2 {r)} / 

(P 1 I 1 + P 2 / 2 ) 2 • (18) 

Here, gi(r) is the usual correlation function of a single two- level system and g^r) that 
of two two-level systems which are dipole-dipole interacting. These correlation functions 
are given in the Appendix. The transition probabilities Pij{r) from period i to period 
j in time r are easily calculated from equation (|9|) and are also given in the Appendix, 
as is Pi, the probability for the occurrence of period i. 

The accuracy of equation (|T8|) is checked for two values of ^3 in figure |2| for a 
coupling constant ReC*3 = — 0.09 A 3 . If the strong laser and the atomic dipole moments 
are perpendicular to the atomic connecting line this corresponds to a local maximum 
at an atomic distance r = 2.7 A, where A is the wave length of the strong 3-1 transition. 
The agreement is excellent and becomes even better if one takes higher orders of C3 into 
account, instead of only the first order as done here. 

The usefulness of equation (P|) is twofold. First, it can be used for a 
fit to experimental data to obtain atomic parameters and mean period lengths. 
Experimentally, correlations functions are in general much easier to determine than 
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Figure 2. g{f) for two dipole-interacting V systems with an atomic distance r = 2.7 A. 
The result of equation (solid line) is compared with a numerical calculation 

(dashed line). The hump after r > 10^3 is clearly visible. The difference for small 
t comes from the the restriction to first-order terms in C3. (a) = 0.3 A3 and 
2 = 0.005 A 3 ; (b) fi 3 = 5.0 A 3 and ft 2 = 0.05 A 3 . 
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period lengths. For such a fit one could also employ numerical solutions for g(r) 
obtained with the quantum regression theorem |34| or with the quantum jump approach 



36, |37], 20]. However, the numerical approach is in general much more sensitive to 
experimental errors than a fit based on an analytic expression. 

The second use of equation (|T8| ) lies in the feasibility to study the behavior of g{r) 
for all parameter values simultaneously, without having to perform many numerical runs 
and possibly overlooking interesting parameter values. This will be demonstrated here 
by exhibiting a possible experimental test of the dipole-dipole interaction. 

4. Possible experimental verification of the dipole-dipole interaction 

For r <C Ti, T 2 , one has (r) = 5^, so that equation ([Tj|) reduces to 

eW = "WW" • (19) 

In particular, one has gi(0) = 0, and (7 2 (0) for angle-averaged detection |38| is given in 



equation ( |A.4|) of the Appendix. Inserting this and the other quantities into equation 



(|P9|) one obtains 



rm = 2P 2 (Al + 2nl)(Al + (ReC 3 y)N 

9{ ) Al(P 1 N + 2P 2 (Al + 2nl)(Al + 2nl + A 3 ReC 3 ))* 1 ) 

with N = (Aj + 2Vll) 2 + A3IC3I 2 + 2A\Re C 3 . To first order in the coupling constant C 3 
this becomes 

1 A 3 (A\ + n\) 2 + n\ „ , . 

By varying Q 3 one has a noticeable change of g{0), which depends on how large ReC3 
is. For small r, the temporal behavior of g(r) is a combination of that of g\ (r) and 
giij) and can be experimentally resolved. Therefore g(0) should be measurable and 
the presence of the dipole-dipole interaction detectable. The deviation of g(0) from 
1/2, i.e. from that for non- interacting atoms, is greatest for small Q 3 and drops off to 
zero for increasing Q 3 . In figure [3] we have plotted g(0) for various values of and 
for ReC 3 = 0.2 A 3 , - 0.1 A 3 , 0AA 3 , - 0.09 A 3 . If the strong laser and the atomic 
dipole moments are perpendicular to the atomic connecting line this corresponds to the 
atomic distance r =1.2 A, 1.7 A, 2.2 A, 2.7 A, respectively. Since ReC 3 oscillatingly drops 
off to zero for increasing atomic distance, its influence on g(0) diminishes for increasing 
atomic distance. As shown in figure |4] the effect is about a factor of two smaller in the 
case of two dipole-dipole interacting two-level systems. 

To show that the temporal behavior of g{r) should be experimentally resolvable we 
have plotted g{j) in figure for various values of Q 3 and for ReC 3 = — 0.09 A 3 , which 
corresponds to an atomic distance of r = 2.7 A for perpendicular laser incidence. 
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Figure 3. g(0) of equation fed) for two dipole-interacting V systems. 
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Figure 4. g(0) of equation (A. 4) for two dipole-interacting two-level systems. 



5. Discussion 

It has been shown for systems with different fluorescent periods that their intensity 
correlation function can be reduced to those of simpler subsystems and to quantities 
that govern the stochastic behavior of the different periods. This is a considerable 
simplification and allows the determination of an algebraic expression for the intensity 
correlation function which usually have to be calculated numerically via the quantum 
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Figure 5. g(j) to first order in C3 for two dipole-interacting V systems with an atomic 
distance of r = 2.7 A and tt 2 = 0.005 A3. 

regression theorem or via the quantum jump approach. 

One of the advantages of an algebraic expression is the feasibility of studying its 
behavior for all parameters simultaneously, without having to make numerical runs 
for different sets of parameters and possibly overlooking interesting values. Another 
advantage is that it is easier to fit experimental data to an algebraic expression and to 
obtain unknown atomic parameters or period durations by such a fit rather than fitting 
to numerically determined expression since the latter procedure is usually much more 
sensitive to experimental and numerical errors. 

It has been shown in this paper that the correlation function for "blinking" systems 
with dark periods shows a "hump" larger than 1 for values of r larger than the correlation 
times of the individual intensity periods. This generalizes a similar result for a V system 
in [0. It has also been pointed out that the existence of extended dark periods may 
considerably enhance the amplitude of Bloch oscillations. 

As an application we have derived an expression for the intensity correlation 
function of two dipole-dipole interacting V systems and have studied its behavior with 
respect to the strength of the dipole-dipole interaction. We have shown that there is a 
significant trace of the dipole-dipole interaction in g{r) for small r which is particularly 
pronounced for small values of the strong driving. This effect of the dipole-dipole 
interaction should be experimentally verifiable for atomic distances of a few wavelength 
of the strong transitions. For two dipole-dipole interacting two-level systems the effect 
is about a factor of two smaller. 
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Appendix A. 



The temporal intensity correlation function for a single two-level system is given by 



1 - e 4 



' 3A > 
cos 7r H sin 7r 

v 4 7 / 



(A.l) 



with 7 = V^ft 2 - A 2 /4, and the mean intensity is 

Aft 2 
1 ~ A 2 + 2ft 2 ' 

We have used the approach of |20| and Eq. (27) of |36| to determine the corresponding 
expression for two dipole-dipole interacting two-level systems and find 



(A.2) 



92{r) = 1 



1 -lAr ( 3A . ^ 
-e 4 C0S7T H sm7T 

2 V 4 7 j 



1 AReC 
'2(A 2 + 2Q 2 )- 



-\At 
e 4 



4ft 2 + 



, , A(A 2 + 2Q 2 )(A 2 



22ft 2 



i6r 



■ r cos 7r 



/512ft 6 + 41A 6 + 2A 2 ft 2 (776ft 2 - 391A 2 ) (A 2 - 6ft 2 )(A 2 + 2ft 2 



64Ar 



47 



to first order in the dipole-dipole coupling constant C. To all orders in C one has 



92(0) 



A 2 + (ReCy 



2A 2 



1 + 



A(A (ImC) 2 -4ft 2 ReC) 
{2Q 2 + A{A + ReC)) 2 



and 



2Aft 2 (2ft 2 + A(A + ReC)) 



z {A 2 + 2ft 2 ) 2 + A 2 Re C{2A + Re C) + A 2 (lm C) 2 
The dipole-dipole coupling constant C is given by |20], [40| 

1 1 



t J sin 7T 
(A.3) 

(A.4) 
(A.5) 



n — ikr 



1 / 

ikr ^ 



cos 2 ■& + 



1 - 3 cos 2 ■& 



(At) 2 i{krf J ^ VJ ' ( A ' 6 ) 

where $ denotes the angle between the dipole moments of the atoms and their connection 
line and k is the wave number of the strong transition. We assume 1? = tc/2 for maximal 
values of C. In Sections [| and £| these expressions are used with A = A 3 , ft = ft 3 and 
C = C 3 . 

For systems with two different light periods and a dark period the relevant transition 
probabilities, P^(t), from period % to period j in time r are found from equation (H). 
For physical reason it is assumed that P02 and P20 vanish. The eigenvalues of the matrix 
B are then /i = and 



Mi,2 = -gCPbi + P10 +P12 + P21) ± 2V^ 01 -P12 - P21) 2 + 4pi pi2 • 



(A.7) 



From this one obtains by means of equation (fLO 



Pn(r) 



P01P21 



M1M2 ^1(^1-^2) 



(?>lo(P21 + Ml) + Pl2(P01 + Ml 
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/^(Ml 



^i(r) 
P22M 



M01M12 



M2) 



'01 



/i 2 , 



/il/i 2 


+ 




P01P21 


+ 


P2ie MlT 


/il/i 2 


A*i(/"i - M2) 


M01M12 


+ 





(P01 + Mi) — 
(Poi + Mi) - 



Pl2Z 



M2T 



yU 2 (/il - M2) 



P2ie 



M2T 



Ml/i2 



M2(Ml 
(Pl2 + P21 + M2) - 



M2) 

P2ie 



(Poi + M2) 
(Poi + M2) 



M2T 



M2(Ml - M2) 



(Pl2 +P21 + Ml) 



and, for r 



00, by means of equation ( |TT 
M10M21 



Po 
Pi 



P10P21 + P01P12 + P01P21 
M01M21 

M10M21 + M01P12 + P01M21 
P01P12 
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(A.8) 
(A.9) 
(A.10) 
(A.ll) 

(A.12) 
(A.13) 
(A. 14) 



M10M21 + M01M12 + M01M21 
For two dipole-dipole interacting V systems with metastable state (figure |l|) the Pij 
have been calculated in f33fl . For this V system the dipole-dipole coupling constant C2 
can be replaced by since A 2 is very small. For simplification we have also put A 2 = 
which does not change the overall results. For zero detuning the transition rates p^ 
between the three periods are then given in |33| as 

2A 3 n 2 2 



Poi 



P10 



P12 



P21 



^3^2 



A 3 



+ ReC 3 



2AI 



2AI 



(Al + 2nj)Ql 



(Al + 2fil)fi§ 



+ ReQ 



(A.15) 
(A.16) 
(A.17) 

(A.18) 



to second order in f2 2 and first order in C3, with the remaining p t j being zero. The 
calculation of the expressions to all orders in C3 and for nonzero detuning is also given 
in p3(] . For the period durations one has T = l/poi> Ti — l/(Pio + M12) and T 2 = l/p 2 i. 
It is noteworthy that Pi becomes independent of Q 2 , up to second order and to all orders 
of C 3 . 
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